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ABSTRACT 

Lopsidedness is common in spiral galaxies. Often, there is no obvious external cause, 
such as an interaction with a nearby galaxy, for such features. Alternatively, the lop- 
sidedness may have an internal cause, such as a dynamical instability. In order to 
explore this idea, we have developed a computer code that searches for self-consistent 
perturbations in razor-thin disc galaxies and performed a thorough mode-analysis of 
a suite of dynamical models for disc galaxies embedded in an inert dark-matter halo 
with varying amounts of rotation and radial anisotropy. 

Models with two equal-mass counter-rotating discs and fully rotating models both 
show growing lopsided modes. For the counter-rotating models, this is the well-known 
counter-rotating instability, becoming weaker as the net rotation increases. The m = 1 
mode of the maximally rotating models, on the other hand, becomes stronger with 
increasing net rotation. This rotating m = 1 mode is reminiscent of the eccentricity 
instability in near-Keplerian discs. 

To unravel the physical origin of these two different m = 1 instabilities, we studied 
the individual stellar orbits in the perturbed potential and found that the presence of 
the perturbation gives rise to a very rich orbital behaviour. In the linear regime, both 
instabilities are supported by aligned loop orbits. In the non-linear regime, other orbit 
families exist that can help support the modes. In terms of density waves, the counter- 
rotating m = 1 mode is due to a purely growing Jeans-type instability. The rotating 
m = 1 mode, on the other hand, grows as a result of the swing amplifier working 
inside the resonance cavity that extends from the disc center out to the radius where 
non-rotating waves are stabilized by the model's outwardly rising Q-profile. 

Key words: instabilities - galaxies: kinematics and dynamics - galaxies: spiral - 
galaxies: structure 



1 INTRODUCTION 

The stellar and/or gaseous discs of spiral galaxies are of- 
ten affected by large-scale asymmetries. About half of all 
late-type galaxies show a lopsided structure that affects the 
whole disc ( |Richter & Sancisi 1994| Haynes et al. 1998). 
Based on near- infrared images of a sample of 149 disc galax- 
Bournaud et al.| ( 2005| find that a large fraction of them 



have asymmetric stellar discs. The strength of the lopsid- 
edness does not correlate with the presence of companions 
but, instead, correlates with the presence of bars and spi- 
ral arms. They explore three different causes for the lop- 
sidedness: galaxy interactions, galaxy mergers, and gas ac- 
cretion. These authors favour the latter explanation, which 
indeed can trigger strong lopsidedness if the gas in-fall is suf- 
ficiently asymmetric. lAngiras et al. (2006) analysed Hi sur- 



face density maps and R-band images of 18 galaxies in the 



* E-mail: Vanessa.Dury@UGent.be 



Eridanus group. All galaxies showed significant lopsidedness 
in their Hi discs. Where the stellar and gaseous discs over- 
lap, their asymmetries are comparably strong. Since the Eri- 
danus group galaxies are more strongly lopsided than field 
galaxies, these authors conclude that tidal interaction in the 
group environment may contribute to generating lopsided- 



ness in disc galaxies. |S waters et al. ( 1999) studied kinematic 
lopsidedness in two spiral galaxies and argued that it may 
be related to lopsidedness in the potential. Alternatively, the 
stellar disc may lie off-centre in the halo's gravitational po- 
tential well and spin in a sense retrograde to its orbit about 
the halo centre (Levine & Sparke 1998). Lopsided structures 
are not the prerogative of disc galaxies alone. In some nu- 
cleated dwarf elliptical galaxies, for instance, the nucleus is 
displaced with respect to the centre of the outer isophotes 
( Binggeli et al.|2000| . While some authors regard these dis- 
placed nuclei as being globular clusters projected close to the 
galaxy photocenter ( Cote et al.||2006| ), in some cases, such 
as the Fornax dwarf elliptical FCC046, there are clear indi- 
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cations that the nucleus is an integral part of the galaxy's 
stellar body and is displaced by a mechanism that affects 
the whole galaxy (De Rijcke & Debattista ||2004 ). 



Sellwood & Valluri ( 1997| used N-body simulations to 



investigate the stability of a family of oblate elliptical galaxy 
models and found that a strong lopsided instability occurs 
in models with small radial anisotropy and strong counter- 
rotation. This instability may cause lopsidedness in stellar 
systems with no or small net rotation, such as dwarf ellipti- 
cals ( De Rijcke fe Debattista|2004| ). |Bertola fc Corsini| | |1999 ) 
compiled a list of several types of counter-rotation in galax- 
ies of different morphological type. They regard counter- 
rotation of stellar vs. stellar discs as the prevailing type of 
counter-rotation since it is the end state of a galaxy with an 
embedded counter-rotating star-forming gas disc. 



Rubin et al. 



( 1992| ) discovered in the SO galaxy 
NGC4550 two distinct stellar components rotating in oppo- 
site directions. Another example is the normal Sab galaxy 



NGC7217 (Merrifield and Kuijken 1994). Vergani et al 



(2007) present the case of NGC5713, a Sbc spiral galaxy in 



which 20 % of the stars are on retrograde orbits. Their data 
suggest that NGC5713 accreted neutral gas from its sur- 
roundings on retrograde orbits. This gas was subsequently 
converted into stars in a counter-rotating disc. Less than 
10 % of all SOs host counter-rotating stellar populations 
while counter-rotating gas is found in roughly one quarter of 
them (Kuijke n et al.|1996| ) . [Kaxmappan Fabrica nt (2001) 
find 4 counter-rotators among a sample of 17 elliptical and 
lenticular galaxies. They find no counter-rotation among 38 
Sa-Sbc galaxies. 

An axisymmetric galactic disc perturbed by a constant 
lopsided halo potential causes a net lopsided distribution in 



the disc, opposite to the perturbation halo potential (Jog 



1997 1999). Tremaine &; Yu (2000) suggest that counter- 



rotation can be produced when a triaxial halo with an ini- 
tially retrograde pattern speed slowly changes to a pro-grade 
pattern speed. This is a situation that probably does not oc- 
cur very often. Moreover, this mechanism requires evolution 
on long timescales of the order of ~ 10 10 yr. Overall, counter- 
rotation is detected only rarely in disc galaxies. Hence, the 
lopsidedness observed in so many disc galaxy is most likely 
not caused by counter-rotation. 

In this paper, we investigate the role played by insta- 
bilities in generating lopsidedness in isolated disc galaxies 
using a semi-analytic matrix method developed by Vauterin] 



& Dejonghe (1996). More specifically, we want to explore 
whether lopsided instabilities can be triggered in fully ro- 
tating disc galaxies. In the next section, we introduce the 
formalism underlying the computer code that we developed 
to analyse the stability of a given dynamical model for a 
disc galaxy. In section [3] we present the unperturbed toy 
galaxy models whose stability is analysed in section^] We 
investigate under what physical circumstances (i.e. degree of 
counter-rotation and orbital anisotropy) lopsided structures 
can spontaneously grow in these disc galaxy models. In sec- 
tion [5] we give a physical explanation for the self-consistent 
growth of lopsided structures, based on the response of in- 
dividual stellar orbits to the growing instability. We sum- 
marise our conclusions in section [6] 



2 SEARCHING FOR INSTABILITIES 

We have developed a computer code to analyse the stabil- 
ity of razor-thin stellar discs embedded in an axisymmetric 
or spherical dark matter halo. The halo is assumed to be 
dynamically too hot to develop any instabilities. This inert 
halo only enters the calculations by its contribution to the 
global gravitational potential. We only consider the stellar 
component of the disc and neglect the dynamical influence 
of gas and dust. 

We describe an instability as the superposition of a 
time-independent axisymmetric equilibrium configuration 
and a perturbation that is sufficiently small to warrant the 
linearisation of the Boltzmann equation. The equilibrium 
configuration is characterised completely by the global po- 
tential Vo(r) and the distribution function fo(E,J), with 
binding energy E and angular momentum J. A general per- 
turbing potential can be expanded in a series of normal 
modes of the form 



V'(r,0,t) = V , (r)e i{md - u,t \ 



(1) 



with a pattern speed R(oo)/m and a growth rate Of (a;), that, 
owing to the linearity of the relevant equations, can be stud- 
ied independently from each other. We write the response 
of the distribution function to a perturbation as: 



f(r,0,v r ,V0,t) = f (E, J) + f'(r,0,v r ,vg,t). 



(2) 



The evolution of the perturbed part of the distribution func- 
tion is calculated using the linearised collision-less Boltz- 
mann equation: 

%-[f',E] = [f y\. (3) 

We rewrite the right-hand side of the last equation as: 
[fo,V,] = -V v /o-V r V" 

dE r d.J 09 y ' 



V r V' • v = 



(5) 



We also know that: 

dV dV 
dt dt " 

With this last identity, equation ([3]) becomes: 

of , a/odv' ., df 0/0 u// fR . 
^T- [f ' E]= 0E ~dT - m dJ )v ■ (6) 

The left-hand side of eq. Q is simply the total time deriva- 
tive of f along an unperturbed orbit. If we integrate eq. 
(|6| along the unperturbed orbits, we immediately obtain 
the response of the distribution function to the perturbing 
potential given by eq. ([!]): 

/'(ro, v ;to) = 



dE 



t 7" / / 1 \ 1 •/ dfo dfo\ f ° T r// \ i(md— cut) j , /^\ 

V (r ;t ) + i(uj — ~ m ~Qj) J V {r)e K j dt.(7) 

The integral in eq. |7| converges if the perturbation disap- 
pears for t — > —00 and is growing sufficiently fast in time 
(3(«)>0). 

Along an unperturbed orbit, the radial coordinate r is 
a periodic function of time with angular frequency u r , just 
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like v r and vq . Because the mean value of vq can be different 
from zero, will be the superposition of a periodic function 
p (t) and a uniform drift velocity cue- 



6 = uj e t + 6 p (t). 



(8) 



We separate the part of the integrand in eq. |7| that is peri- 
odic with frequency uo r from the aperiodic part and expand 
it in a Fourier series: 



i(mujQ—uj)t \ ^ t ilijj r t 



The coefficients Ii are given by 



e at, 



(9) 



(10) 



where the integration extends over half a radial period, start- 
ing at apocentre at t = 0. 

The first term of the right-hand side of equation |7]) is 
rewritten as: 



J — oo 

(11) 



' t ri t , \ QjO i(m0 

dE V{r °' t0) = dE 6 



dfo o i{md Q -ujt Q ) 

dE° 



The orbits are integrated in the unperturbed potential 
using a leapfrog method. Instead of using E and J, orbits 
in the unperturbed potential are catalogued by their apoc- 
entre and pericentre distances, denoted by r+ and r_, re- 
spectively. The sense of rotation is indicated by the sign of 
r_. The grid of the orbit catalogue is given by (r+,r_) £ 
[0,r max ] x [— r+,r+]. We use a grid of 100 x 200 cells. For 
every orbit, cj r and uq are determined, the Fourier expan- 
sion (|9| is performed up to the order of / max = 50, and the 
coefficients are stored. If we want to calculate the response 
of the distribution function in the point (i*o,vo) in phase 
space at time to , we choose an orbit from the catalogue with 
the correct integrals of motion (i.e. r+ and r_) but passing 
through its apocentre at t = so the actual orbit has an 
offset in time t(ro) and azimuth p (ro) that must be taken 
into account. Therefore we also store a tabulation of t(r) 
and p (r). 

The response of the distribution function to the pertur- 
bation now assumes the following concise form: 



/'(r , v ;to) 



i(m0 o -ut o ) 



i= 



luj r + mujQ — UJ 



Because we are searching for instabilities, and thus £y(u;) > 
0, we do not have to be concerned by the presence of reso- 
nances. We can compute the perturbed density by integrat- 
ing the perturbed distribution function over the velocities 
up to the escape velocity. For self-consistent perturbations, 
the gravitational potential produced by this response den- 
sity should equal the original perturbing potential. Using the 



matrix method developed by Vauterin & Dejonghe (1996), 
the search for a self-consistent mode of order m is reduced 
to an eigenvalue problem. We adopt the same basis set of 



100 




Figure 1. Rotation curve of the unperturbed potential. The dot- 
ted line represents the contribution of the disc. 



The response density generated by each of these basic poten- 
tials can be expanded in terms of the basic density distribu- 
tions. This gives rise to a matrix C(uj) that contains the co- 
efficients of these expansions. Self-consistent perturbations 
have the unique property that they have a pattern speed and 
growth rate uj for which C{uS) has an eigenvalue A = 1. The 
corresponding eigenvector contains the coefficients of the ex- 
pansion of the response density in terms of the base set of 
density distributions. Thus, for a self-consistent perturba- 
tion of order m we are left with the numerical search within 
the complex plane for a value of uj for which C(u>) has a 
unity eigenvalue, something that can be accomplished very 
efficiently using bisection. For a more detailed description 
of the method, we refer the reader to Vauterin & Dejonghe 
([1996}. 



3 THE UNPERTURBED MODELS 

A dynamical model for a disc galaxy embedded in a dark 
halo is specified by the global potential of the system and the 
distribution function of the stellar disc. For all the models 
that will be analysed in this paper, we have chosenb^ the 
same unperturbed potential and unperturbed density but 
different distribution functions. The unperturbed potential 
in the plane of the disc is given by: 

GM GM 



V (r) = 



+ 



(13) 



with M = 1O 1O M0 and r expressed in kpc. This potential 
produces a rotation curve that rises near the centre and be- 
comes flat further out in the disc, as can be seen in Fig. [l] 
For the unperturbed density profile we choose an approxi- 
mately exponential profile: 



po — ae 



-1.3^0.2+r 2 



(14) 



potential-density couples as Vauterin & Dejonghe (1996) 



The mass of the disc, and thus the mass of the halo, is 
determined by a. We impose an outer limit at r ma x = 6 
kpc. We have chosen a so that the proportion H/D of the 
total mass inside the outer radius r max for the halo and 
the disc is about 2.5. The disc is truncated by demanding 
that the distribution function fo(E, J) is zero for orbits that 
venture outside this outer limit. The distribution function is 
written as a linear combination of a basic set of distribution 
functions. The coefficients of this expansion are determined 
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Figure 2. Unperturbed distribution functions (left) and unperturbed moments (right) of a few analysed models: (a) /q'°, (b) /q' 25 '°, 
(c) /q' 5 '°, (d) /q -5 ' 0-5 , (e) /q -5 ' 1 , (f) /q' 1 - For the distribution functions the density is decreasing linear from black to white. 
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by a least-square fit to the density profile given by eq. ( 14 ) 



The dominant m — 1 mode we find in strongly counter- 
rotating models is obviously the well-known counter- 
rotation instability. This instability has been known since 



Zang &; Hohl| (|1978|) and has been stud ied analytically 



(Sawamura 1988 [Palmer &; Papaloizou 1990 Tremaine 
2005) and using N-bod y simulations (|Merritt &; Stia velli 



1990l|Levison et ed.|1990||Sellwood fe Merritt |1994||Sellwood 



|fe Valluri||1997| ). jSellwood fe Valluri| (1997) donot detect 
it in systems rounder than E6, which is mainly the result 
of their rounder models being stabilised by a higher radial 



In all the models, the error on the fit to the mass density 
never exceeds 1% of the central value. The contribution of 
the disc to the rotation curve is shown in Fig. [I] 

We start with a strongly tangentially anisotropic model 
that is know n to develop strong spiral arms ( |Vauterin &; De- 
jonghe|l996| . Its distribution function f^°{E, J) is shown in 
the left column of Fig. [2^l. Most stars in this model populate 
nearly circular orbits that rotate in one direction. The un- 
perturbed moments are shown in the right column of Fig.[2^L. 

Using the distribution function fQ ,0 (E, J) = /q ,0 (E, — J), pressure. On the other hand, Merritt &; Stiavelli (1990) find 
which has the same unperturbed density as /q ,0 (E, J), one 
can easily construct counter-rotating discs that all have the 
same density distribution. In order to place a fraction x of 
all stars on retrograde orbits, the following distribution func- 
tion can be used : 



f^(E, J) = (1 - x)f ' u (E, J) + xft u {E, J) 



(15) 



In Fig. [2J3, we plotted the distribution function and unper- 
turbed moments of a model in which 25 per cent of the 
stars counter-rotate. In Fig. |2|;, we present the properties of 
a model consisting of two equally massive counter-rotating 
discs. 

In the left column of Fig. [2^, we show the distribution 
function / °- 5,1 (£,J) of a very radially anisotropic model 
that is known to produce a strong bar instability ( Vauterin fe| 
Dejonghe 1996). Using this model we can construct models 
with varying degrees of radial anisotropy, as follows: 



f°^(E, J) = (1 - y)/r' u (E, J) + yfr'\E, J). 



0.5,0 / 



rO.5,1/ 



(16) 



In Fig. [2jl, we plotted the properties of a model consisting 
of an even mix of /q 5 ° an d /o" 5 ' 1 (i-e. y = 0.5). We can also 
construct a similar family without counter-rotating discs: 



ff v (E, J) = (1 - y)f '"(E, J) + yft L (E, J). 
In Fig. [2J, we plotted the distribution function of /q 



(17) 



•0,1 



4 LOPSIDED INSTABILITIES 

4.1 The influence of counter-rotation 

In Fig. [3] we show the growth rate and pattern speed of the 
fastest growing m — 1 and m = 2 modes as the fraction 
x of counter-rotating stars for the models of the /q '° fam- 
ily changes. Clearly, the fully counter-rotating /q' 5 ' model 
has a dominant non-rotating m — 1 instability whereas the 
m — 2 instability is virtually absent (actually, two mirror- 
image m = 2 modes with the same growth rates and oppo- 
site pattern speeds occur in the counter-rotating model). As 
the degree of counter-rotation diminishes, the growth rate 
of the m = 1 instability slowly declines while it picks up a 
small pattern speed. At the same time, the m — 2 instability 
increases in strength. As long as more than one quarter of 
all stars are on retrograde orbits, the model is dominated 
by a slowly rotating lopsided mode. Unexpectedly, around 
x ~ 0.08, the nature of the m — 1 mode changes. Another, 
this time rapidly rotating, lopsided instability becomes the 
fastest growing m—\ mode. Its strength increases together 
with that of the m — 2 mode as the degree of counter- 
rotation vanishes. 



lopsidedness developing in systems as round as El but with 
negligible radial pressure. Partial rotation only introduces a 
pattern speed in an otherwise purely growing instability. 

The counter-rotating bars we found, were also re- 
port ed by jSellwood fe Merritt | fl!994|), |Levison et al.| ( |1990 ) 
and |Friedli| fl!996| ). |Sellwood & Valhiri| ( |1997| ) argued that 
they were the result of non-linear orbit trapping in finite- 
amplitude spiral disturbances. The fact that we found them 
shows that they are formed through linear instabilities. 

Recently, some authors found a lopsided instability in a 



normal differentially rotating galactic disc (iSaha, Combes & 



Jog 2007). The lopsided pattern precesses in the disc with 
a very slow pattern speed with no preferred sense of pre- 
cession. The weaker m = 1 mode that we found in the 
rotating model has a certain sense of rotation and bears 
strong resemblance to the so-called eccentricity instability 
that occurs in gaseous and stellar near-Keplerian discs or- 
biting a central massive object ([Adams, Ruden &; Shu| 1989 



|Shu et al.|1990l|Noh, Vishniac Cochran|1991| |Taga &; lye 
1998| |Lovelace et al.|1999|j Jacobs ■ &: Sellwood|2001||Baconet 
|al.|2001||Salow fe S t at ler||2 004 1 ) . | Lovelace et al.| ( |1999D con- 
structed models for disc galaxies with exponentially declin- 
ing surface density profiles embedded within a spherically 
symmetric dark halo. These authors found the inner regions 
of such systems rapidly develop a trailing one-armed spiral 
wave, even if the mass of the central object is small. The 
first N-body example of a rotating lopsided instability was 
found by Sellwood (1985]) in a mass model of our Galaxy 
without a halo component. | Evans &; Read] ( |1998"t examined 
the global stability of stellar power-law discs. They report a 
similar rotating lopsided pattern in cut-out power-law discs, 
but found no growing non-axisymmetric modes in the fully 
self-consistent power-law discs. We provide the first theoret- 
ical evidence, based on a thorough mode-analysis of a suite 
of self- consistent dynamical models for disc galaxies embed- 
ded in a dark halo, that the eccentricity instability can also 
occur in the fully prograde stellar discs of spiral galaxies 
without an additional massive central component, such as 
a compact bulge or super-massive black hole, and without 
introducing an unresponsive central region in the disc. 

4.2 Anisotropy 

We now study the behaviour of the self-consistent m — 1 and 
m — 2 as the fraction y of stars on radial orbits changes. The 
variation of the growth rate and pattern speed is presented 
in Figs. [4^ and jij) for the models /o -5 '^ (black) and /q ,v 
(red). We start from model /q -5 ' that is known to develop a 
strong non-rotating m — 1 instability and two twin counter- 
rotating bar instabilities. As y, and thus the fraction of stars 
on radial orbits, increases, the m — 1 and m — 2 instabilities 
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Figure 3. Growth rate (a) and pattern speed (b) for the strongest m = 1 (solid line) and m = 2 (dotted line) instabilities as the 
fraction x of counter-rotating stars changes. Discontinuities appear when the nature of the instability changes. The perturbed density 
p'{p = po + ep f ) and velocity field v'(v = vo + ev') is plotted out to a radius of 2.5 kpc: (c) rotating one-armed spiral for /q'° , (d) rotating 
two- armed spiral for /q'°, (e) non-rotating lopsided mode for /q" 5 '°, (f) one of a pair of counter-rotating bars for /q" 5 '°. Over-densities 
are coloured white- red, under-densities black-blue. In panels (c) and (e), we also indicate the radius within which a growing m = 1 mode 
is expected based on density wave theory (see section [5] for a detailed discussion). 




Figure 4. Growth rate (a) and pattern speed (b) for the strongest m = 1 (solid line) and m = 2 (dotted line) instabilities as the fraction 
y of stars on radial orbits changes for the non-rotating (black) and rotating models (red). The perturbed density p' (p = po + ep') and 
velocity field v x (v = vq + ev x ) is shown for model /q" 5 ' 1 in (c). Over-densities are coloured white-red, under-densities black-blue. 



rapidly stabilise. The lopsided mode is the first to stabilise, 
at y ~ 0.05. The model becomes fully stable against both 
m — \ and m = 2 modes at y ^ 0.2. For y > 0.6, the 
model develops a non-rotating bar instability that becomes 
stronger as radial anisotropy increases. The perturbed den- 
sity and velocity field of the bar of the /q' 5 ' 1 are shown in 
Fig. 13;. If we start from model /q'°, the m = 1 and m = 2 



instabilities also stabilise with increasing anisotropy. The 
m — 2 instability, is the last one to stabilise and the model 
becomes stable for y > 0.6. 

From this exercise, it is clear that the mechanism that is 
responsible for triggering the lopsided mode in the counter- 
rotating model relies heavily on virtually all stars moving 
on near-circular orbits. Even a relatively small contribution 
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of stars on radial orbits makes it impossible for the disc to 
develop the m—\ mode. 

4.3 Perturbed line-of-sight velocity fields 

To allow for a direct comparison of our models with ob- 
servations, we present in Fig. [5] the perturbed line-of-sight 
velocity fields for the extreme cases of exact counter-rotation 
(top row) and full prograde rotation (bottom row) . The first- 
order perturbation to the velocity v', corresponding to the 
perturbation f to the distribution function /o, is given by 



/ v/'dv - p'vo 



Po 



(18) 



This velocity perturbation is then projected onto the sky in 
order to obtain the perturbation on the line-of-sight veloc- 
ity denoted by v[ os . The velocity fields in Fig. [5] have an 
inclination i = 50° and a viewing angle <j) ranging from 0° 
to 90°. The velocity fields are only plotted out to a radius 
of 2.5 kpc since only inside this region is the density pertur- 
bation noticeable. Schoenmakers et al.l (|1997| showed that 



if the potential contains a perturbation of harmonic number 
m then the line-of-sight velocity field contains n = m — 1 
and n — m + 1 terms. Thus as expected, we can see n — 
and n = 2 components in the residual velocity fields but 
there is no strong direct resemblance with the two galaxies 



discussed by S waters et al.| ( |1999| , one of which shows a dom- 
inant n — 2 velocity perturbation while the other exhibits 
predominantly a n = structure. Of course, the residual 
velocity fields plotted in Fig. [5] are calculated using linear 
perturbation theory. Instabilities in real galaxies are likely 
to be in the non-linear regime. Therefore, this comparision 
is intended to be indicative not definitive. 



5 PHYSICAL INTERPRETATION 

In order to unravel the physics behind the two distinct m = 1 
modes found in section [4] we will study the orbits of stars 
that move in the global perturbed potential: 



V(r,e,t) = V (r) + eV'(r)e i(m6 - 



Ot) 



(19) 



In order to simplify the interpretation, we keep the ampli- 
tude of the perturbation fixed, i.e. we set ^s(ui) — 0, and 
only consider its pattern speed, R(oo)/m. The prefactor e is 
determined by requiring that the maximum difference be- 
tween the perturbed and the unperturbed density nowhere 
exceeds 10 % of the unperturbed density. We then numeri- 
cally evolve an ensemble of stars in the perturbed potential 
using a leapfrog integrator. The goal is to see how the orbits 
are affected by the perturbation and which perturbed orbits 
help support the perturbation. For brevity, we will hence- 
forth refer to the m = 1 instability of the /q -5 ' model as 
the "non-rotating lopsided mode", and to that of the /q'° 
model as the "rotating lopsided mode" . 



5.1 Non-rotating lopsided mode 

In Fig. [6] we present some perturbed stellar orbits that oc- 
cupy the region in which the density perturbation is maxi- 
mal. They all share the same energy but are characterised by 
a different angular momentum. We can easily recognise two 





(a)" 



■(b)- 



Figure 6. Stellar orbits in the potential well of the non-rotating 
lopsided mode. All orbits have the same energy but different an- 
gular momentum: Angular momentum increases in figure (a) from 
the black orbit (J = 0) to the cyan orbit (J = J c ) and in figure 
(b) from the cyan orbit (J = J c ) to the green orbit (J = 1.9 J c ). 




Figure 7. (a) Butterfly orbit, (b) Loop orbit. To see the evolution 
in time, we indicated some parts of the orbit in colour. The first 
part red, the second part green and the last part blue. 



general orbit families : butterflies and loops. In Fig. [T^l we 
show a more detailed picture of a butterfly orbit. A butter- 
fly can be viewed as a librating elliptical orbit with variable 
eccentricity. At the turning points of the libration, the orbit 
ellipticity becomes zero (this is evidenced by the red, green, 
and blue ellipses in Fig. |7^). As the angular momentum of 
the orbit is increased up to a critical point J c , the two in- 
stances of zero ellipse orbit ellipticity coincide and the orbit 
fills an elliptical region. For still higher angular momentum, 
libration becomes rotation and the orbit becomes a loop, as 
can be seen in Fig. ^6jp & Fig. [Ffc. Remar kably, both orbit 
families are also found by Jalali & Ranee (2000) for a disc 
galaxy model with a lopsided potential that is of Stackel 
form in elliptic coordinates and with two separate strong 
density cusps. They also found two other orbit families, nu- 
cleophilic bananas and horseshoe orbits. It is clear from their 
formulation that these two orbit families are associated with 
the cusps having diverging central densities which is why we 
do not find them in our models. 



5.1.1 Linear regime 

The mechanism that causes the non-rotating m — 1 mode 
is now clear : an infinitesimal m — 1 perturbation will cause 
near-circular orbits to become somewhat more elliptic and to 
shift towards the slight overdensity, thus adding to this over- 
density, which in turn causes other orbits to become more 
elongated and to shift, and so on. This was already evident 
from section 3.3 3(a) of Binney & Tremaine ( |1987| who used 
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Figure 5. Residual line-of-sight velocity field v[ (with i;i os = i>i os ,o + ev l os ) f° r ^ ne m — 1 instabilities of (a) / " ' and (b) / ' . All 
fields have an inclination of 50° and a viewing angle 4> varying from 0° to 90° in steps of 30°. Positive velocities are coloured white-red, 
negative velocities black-blue. 
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Figure 8. (a) The sign of the quantity C2(r) (solid line) roughly traces that of the density perturbation (dotted line). Perturbed orbits 
(green) are displaced from the unperturbed orbit (red orbit): orbits with C2 > are shifted into the direction of the overdensity (b), 
those with C2 < are shifted into the opposite direction (c). 



perturbation theory in combination with the epicyclic ap- 
proximation to calculate the response of near-circular orbits 
to a general m-armed perturbation. Since this perturbation 
has no Lindblad resonances, the condition that the sign of 
the quantity Cz(r) = dV'(r)/dr + 2V'(r)/r (represented by 
the red curve in Fig. [8^l) traces that of the orbit displace- 
ment is fulfilled everywhere within the stellar disc (see eq. 



(3-120b) of |Binney fe Tremaine| fl!987| )) and the perturbed 
near-circular orbits will all strengthen the m = 1 mode. Tan- 
gential orbits with C2 > are shifted into the direction of 
the overdensity (Fig.j8j)), those with C2 < are shifted into 
the opposite direction (Fig. |8j:). As a consequence, the sign 
C2 also roughly follows that of the density perturbation. 

With a maximum density contrast of 10%, the orbits 
plotted in Figs. [6] and [7] are not per se in the linear regime. 
This was done, since we here only wish to illustrate the 
mechanism that induces the instability to grow, for clar- 
ity : a truly infinitesimal perturbation would result in an 



equally infinitesimal and therefore nearly invisible shift of 
the tangential orbits. 

It is instructive to interpret this instability not only 
on the level of stellar orbits but also in terms of density 
waves. As shown in Palmer's book on dynamical insta- 



bilities ^Palmer 1994), a razor-thin disc consisting of two 



equal-mass counter-rotating stellar populations may develop 
purely growing one-armed (m — 1) WKB waves if the 



Toomre parameter, Q — 
stability criterion 



3.36Gp ' 



K(r) 2 - fi(r) 



2 * 



fails to satisfy the local 



(20) 



This is very similar to the well-known criterion Q ^ 1 for 
local stability against purely growing m — waves. The pre- 
cise form of this stability criterion depends on three approx- 
imations, (i) the WKB approximation for tightly wound spi- 
ral waves, (ii) the epicycle description for near-circular stel- 
lar orbits, and (Hi) a Gaussian distribution function. None 
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Figure 9. In a counter- rotating disc, purely growing WKB waves 
may develop when Q(r) (dashed line) ^ -, ff n , , 2 (solid line). 

K{r) ^H r j ' 



The quantity 



^( r ) / 2_0( r )2 is plotted using a full line. The Q— 
profiles of several counter-rotating models, with increasing radial 
anisotropy as quantified by the parameter y in eq. ( |16| , are over- 
plotted (see legend in the figure). Models with y > 0.06 are stable, 
both according to Palmer's stability criterion and our mode anal- 
ysis. 



of these necessarily applies to the models presented in this 
paper. Nonetheless, we used the Q(r), ft(r), and Q(r) pro- 
files of the fully counter-rotating /q ' 5,y models to evaluate 
this criterion as a function of radius, where Q depends on 
the parameter y via the radial velocity dispersion (Fig. |9j). 
We only check the fully counter-rotating models since only 
in this case can an analytical stability criterion be derived. 
It is clear from Figs. [3|e) & [9] that the region where, ac- 
cording to the local stability criterion, purely growing waves 
may develop roughly coincides with the region in which 
the non-rotating lopsided mode resides (i.e. inside a radius 
r < 2.1 kpc). 

As a further test, we checked whether the model that, 
according to our mode-analysis, is the first to be stabilized 
by increasing radial anistropy, is also stable according to 
Palmer's criterion. This is done by increasing the parameter 
y in eq. (16). It is clear from Fig.|9]that the Q(r) profile of 
the first stable model according to our mode-analysis, the 
one with y « 0.06 (see fig. |4|l), is also very close to the 
line of stability according to Palmer's criterion. For y ^ 0.1, 
the models are definitely stable both according to Palmer's 
criterion and our mode-analysis. Note that Palmer's stability 
criterion is a sufficient one : satisfying it implies stability, 
not satisfying it does not necessarily imply instability. Given 
the reasonable agreement between Palmer's analysis and our 
mode analysis concerning the line of stability and the spatial 
extent of the m = 1 pattern, we are led to the conclusion 
that the purely growing lopsided mode we find in counter- 
rotating discs is caused by a local Jeans-type instability. 



5.1.2 Non-linear regime 

Figs. [6] and [7] re veal a feature that is not captured by our lin- 
ear mode-analysis. Once the perturbation is strong enough, 
radial orbits are deformed into the new family of butter- 
fly orbits whose centres of gravity are also shifted in the 
direction of the over-density. Thus, this orbit family may 




Figure 10. The change in the density distribution of an ensemble 
of orbits, 5p, as defined by eq. ( |21| >, due to the non-rotating m = 1 
mode. Panel (a) : an ensemble of circular unperturbed orbits; 
panel (b) : an ensemble of radial unperturbed orbits. Negative val- 
ues of Sp are coloured black/blue; positive values in white/yellow. 
The contours of the perturbed density of the m = 1 mode are 
also plotted (solid line: over-density, dotted line: under-density). 
See text for a discussion of this figure. 



potentially contribute to the m — 1 perturbation but only 
after the amplitude of the perturbation has become large 



enough, see also Jalali & Ranee (2000). For an infinitesimal 



perturbation, radial orbits do not contribute to the growth 
of the instability (their presence is even detrimental to the 
instability's growth in the linear regime, see below). 

Further evidence for this mechanism can be gleaned 
from the following exercise. A test particle with initial con- 
ditions (r , vo) on an orbit in the unperturbed system with 
potential Vo(r) over time fills a certain area, which, in gen- 
eral will have the form of an annulus (with circular and 
straight-line orbits as extremes). For the same initial con- 
ditions (i*o,vo), the test particle's orbit in the perturbed 
potential V(r) will fill a differently shaped area. The differ- 
ence between the density distribution of the perturbed orbit 
and that of the unperturbed orbit gives an idea of how the 
density of a galaxy made up of an ensemble of unperturbed 
orbits will change under the influence of the given pertur- 
bation. In Fig. [lO] we show how the density distribution of 
an ensemble of circular orbits (lOi) or radial orbits (10d) 
changes due to the perturbation given by eq. (19); i.e., we 
plot the quantity 



Sp(r,6) = (Porbit,pert(r, 6) - Porbit,unpert(r, 6)) 



(21) 



where the sum runs over an ensemble of circular unper- 
turbed orbits with different radii and with the phases of 
the starting points of the orbit integrations distributed uni- 



formly over the interval [0,27r[ (circles in Fig. 10 1 and spokes 
in Fig. 10 3 are caused by the finite number of orbits). Neg- 
ative values of 5p are coloured black/blue; positive values 
in white/yellow. Clearly, the regions where 5p > coincide 
with the over-densities of the m — 1 mode (full line con- 
tours) whereas the regions where 5p < coincide with the 
mode's under-densities (dotted line contours). It is obvious 
from Fig. |10| 3 that the radial orbits are not nearly as co- 
operative despite the fact that they are slightly displaced 
towards the inner overdensity. As an ensemble, they do not 
react to the imposed perturbation in a way that would tend 
to strengthen it since their long axes are oriented perpendic- 
ularly to the inner lopsidedness leading rather to a m = 2 
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feature. The (near-) circular orbits are clearly the backbone 
of this m = 1 mode and, as we know from Fig. [4^,, even a 
smidgen of stars on radial orbits is enough to stabilise the 
system against this instability. 



5.2 Rotating lopsided mode 

In Fig. [TT] we present five perturbed stellar orbits that oc- 
cupy the region in which the density perturbation due to 
the rotating lopsided mode is maximal. The orbits are plot- 
ted in a reference frame that rotates with the pattern speed 
of the m — 1 mode. Again we observe loop orbits that are 
displaced into the direction of the over-density and that can 
support the lopsided structure (Fig. &JTTJ)). At larger 
radii, the stellar orbits attain a banana shape in this refer- 
ence frame. These orbits occupy a region outside the mode's 
main under-density and seem to be connected to the one- 
armed spiral (Fig. [TTJ:, [TTJi fclTT^). Unlike in the previous 
case, where aligned loop orbits were almost solely respon- 
sible for creating the instability, here the two orbit families 
fulfil different tasks. The aligned loop orbits support the in- 
ner lobes of the lopsided structure whereas the banana orbits 
make up the outer one-armed spiral. 



5.2.1 Linear regime 

For the lopsided mode with pattern speed Q p , the quantity 
C2(r) from Binney & Tremaine (1987) becomes: 



C a (r) 



- m 2 (Q — Q p 



dV'jr) flV'(r) 
dr r(fi-fip) 



• (22) 



The sign of Ci again traces that of the orbit displacement 
and, with less fidelity, of the density perturbation (Fig. 
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Orbits inside the corotation radius are shifted into the direc- 
tion of the overdensity (Ci > 0) (Fig. 



12 



d). At corotation, 



C2 diverges. Orbits near corotation become banana orbits. 
The corotation radius of the mode is at r c = 0.8 kpc, which 
coincides with the position of the one-armed spiral. Outside 
corotation, the sign of C2 changes and orbits are shifted into 
the other direction (Fig. 12:). 



The interpretation of this instability in terms of waves 
is somewhat subtle. In the case of a bar instability, the radial 
extent of the m — 2 pattern is determined by the largest ra- 
dius out to which the most slowly rotating wave that avoids 
having an inner Lindblad resonance (ILR) can travel. This 
constraint sets the size of the resonance cavity within which 
the pattern can grow by swing amplification. However, one- 
armed waves do not have an ILR, irrespective of their pat- 
tern speed. The radial extent of the m—\ pattern is set by 
the largest radius out to which non-rotating wave packets 
can travel. This radius is determined by the disc's disper- 
sion relation. We used the WKB dispersion relation (eq. (6- 
46) of|Binney fe Tremaine| fll987| or eq. (12.82) of |Palmer 
(19941) together with the Q(r), ft(r), and Q(r) profiles of 
the / '° model to estimate the region in which non-rotating 
one-armed waves are allowed to exist. The Q(r) profile of 
the fully rotating /q '° model rises outwardly, limiting the ex- 
tent of non-rotating waves to some finite radius, well within 
the disc. In the case of the /q'° model, for a zero pattern 
speed, the dispersion relation has two branches (the long- 
wave and the short-wave branch) of incoming and outgoing 




Figure 13. The change in the density distribution due to the 
rotating m = 1 mode. Panel (a) : an ensemble of circular un- 
perturbed orbits; panel (b) : an ensemble of radial unperturbed 
orbits. Negative values are coloured black/blue; positive values 
in white/yellow. The contours of the perturbed density of the 
m = 1 mode are also plotted (solid line: over-density, dotted line: 
under-density) 

[width=17cm] 

leading and trailing waves for radii smaller than approxi- 
mately 1.7 kpc. This sets the dimension of the resonance 
cavity within which the pattern can grow through swing am- 
plification. All m = 1 waves can propagate into the galaxy 
center where incoming trailing wave packets are reflected as 
outgoing leading wave packets, closing the feedback loop. 
Moreover, a more steeply rising rotation curve will suppress 
the m = 2 mode while leaving the m = 1 mode, which has 
no inner Lindblad resonance, largely intact. 

Thus, following Evans &; Read ( 1998| ) who already pro- 
posed this mechanism as the cause of the m = 1 modes 
found in cut-out power-law discs, we propose swing amplifi- 
cation as the physical cause of the one-armed mode in this 
rotating model. 



5.2.2 Non-linear regime 

In Fig. [13] we show how the density distribution, i.e. the 
quantity Sp, defined by eq. (21), of an ensemble of circular 
(panel (a)) and of radial (panel (b)) orbits changes due to the 
rotating m — 1 mode. We now see that the circular orbits 
and the radial orbits both support the lopsided structure. 
This is because radial orbits as well as the inner circular 
orbits become loop orbits that are aligned with the lopsided 
mode. If we change the fraction of stars on radial orbits, we 
have a slower stabilisation of the m — 1 mode than in the 
non-rotating models (see Fig.[4^). 

We have integrated a number of stellar orbits that all 
start at the corotation radius r — r c with phases = 
45°, 90°, ...225°. The results are shown in Fig. [HI where we 
also showed the effective potential 

V eB (r,0) = V(r,0)- i(9(w))V. 



(23) 



All banana orbits evolve around the stable Lagrange point 
of the system which is located at the local maximum of the 
effective potential at a phase angle 0l- In the corotating 
reference frame, a star at position (r c , 6l) will remain there 
forever. Orbits started at corotation but with different phase 
angles will trace out a banana-shaped curve. The further 
away from the equilibrium point (r c , 0l) an orbit is started, 
the more elongated the banana becomes. The orbit started 
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Figure 11. Stellar orbits in the force field of the rotating lopsided mode, as seen in a reference frame rotating with the same pattern 
speed as the lopsided mode. We find two orbit families : (a) & (b) loop orbits, and (c), (d) & (e) banana orbits. 



(a) 




(b) 




(c) 




Figure 12. (a) The sign of the quantity C2(r) (solid line) roughly traces that of the density perturbation (dotted line). Perturbed orbits 
(green) are displaced from the unperturbed orbit (red orbit). Orbits with C2 > are shifted into the direction of the overdensity (b) and 
those with C2 < are shifted into the opposite direction (c). 




Figure 14. A few banana orbits in a frame rotating with the 
lopsided mode. All orbits start from co-rotation with the same 
energy and angular momentum, but at a different angle 0: white 
(0 = 45°), purple (0 = 90°), cyan (0 = 135°), yellow (0 = 180°), 
black (0 = 225°). The effective potential is also shown (increasing 
from white to blue). 



at (r c , 180° is circular in the corotating reference frame. 
The origin of the one-armed spiral now becomes more clear. 
A weak, rotating lopsided perturbation will catch stars close 
to the (r c , 0l) Lagrange point on banana orbits. These then 
spend more time near the stable Lagrange point, creating 
an over-density there, and less time near the diametrically 
opposite point, creating an under-density there. Thus, the 
banana orbits help strengthen the lopsided mode, causing 
more stars to be trapped in banana orbits, and so on. 



6 CONCLUSIONS 

Using a toy dynamical model, we have investigated the prop- 
erties and causes of dynamical instabilities that may cause 
lopsidedness in disc galaxies. 

We found the well-known counter-rotation instability to 
be the dominant mode in disc galaxy models with strong 
counter-rotation and small radial anisotropy. It is much 
stronger than any m = 2 mode we found in these mod- 
els. The strength of this mode diminishes in models with 
less counter-rotation and it eventually disappears in fully 
rotating models. These, however, develop a different type of 
lopsided mode, that becomes stronger as rotation increases, 
although it is always much weaker than the m — 2 spiral-arm 
mode. This instability bears strong resemblance to the ec- 
centricity instability that is known to occur in gaseous and 
stellar near-Keplerian discs orbiting a central massive ob- 
ject. We provide the first theoretical evidence, based on a 
thorough mode-analysis of a suite of dynamical models for 
disc galaxies embedded in a dark halo, that the eccentricity 
instability can also occur in the fully pro-grade stellar discs 
of spiral galaxies. 

By integrating the orbits of an ensemble of stars in the 
perturbed potential of the two extreme cases of full counter- 
rotation on the one hand and full prograde rotation on 
the other hand, we investigated the physics underlying the 
counter-rotation and eccentricity instability. The counter- 
rotation instability grows by changing near-circular orbits 
into aligned loop orbits that help maintain a lopsided struc- 
ture. If the non-linear regime, radial orbits are changed into 
butterfly orbits. In the case of the eccentricity mode, both 
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radial and tangential orbits become aligned loops in a coro- 
tating reference frame. In the non-linear regime, orbits near 
corotation are trapped into resonance and describe banana- 
shaped figures in a corotating frame. They help support the 
characteristic one-armed spiral. 

In terms of density waves, the counter-rotating m = 1 
mode is most likely due to a purely growing Jeans-type in- 
stability. An approximative analytical criterion for local sta- 
bility, akin to Toomre's stability criterion for axisymmet- 
ric waves, can be employed to estimate the region in which 
purely growing one-armed waves may develop. This estimate 
roughly coincides with the region in which the non-rotating 
lopsided mode is observed to reside. The rotating m — 1 
mode, on the other hand, grows as a result of the swing 
amplifier working inside the resonance cavity that extends 
from the disc center out to the radius where non-rotating ro- 
tating waves are stabilized by the model's outwardly rising 
Q-profile. Rotating waves are confined to even smaller radii 
so the non-rotating waves effectively set the outer boundary 
of the resonance cavity. 

Many disc galaxies show a noticeable m = 1 perturba- 
tion besides the dominant m — 2 spiral-arm mode but only 
very few of them show any counter-rotation. The rotating 
lopsided mode we identified in fully prograde disc models 
therefore forms an attractive explanation for this observed 
phenomenon. 
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